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Abstract 

A symbolic computation technique is developed to calculate adiabatic evolution 
equations for parameters of the perturbed DNLS/MNLS solitons obtained by the 
recently developed direct perturbation theory [X.-J. Chen and J. Yang, Phys. Rev. 
E 65, 066608(2002)]. Effects of the intrapulse Raman scattering, third-order group 
velocity dispersion, and narrow-banded filters with nonlinear gain on the MNLS 
solitons are studied as examples. 
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1 Introduction 



It has been well known that exactly integrable nonlinear differential equations 
support soliton solutions which travel stationarily and collide elastically. It was 
also known that the physical situations giving rise to exactly solvable equations 
are highly idealized. Perturbations violating their integr abilities actually exist. 
If these perturbations are small enough, their influence on solitons can still 
be known analytically by perturbation theories for solitons(see, e.g., [1,2,3]). 
In the picture depicted by perturbation theories, the lowest approximation is 
an adiabatic solution in which the soliton keep its profile unchanged while its 
parameters such as amplitude, velocity, initial center and initial phase may 
evolve adiabatically[l,2]. Perturbations can also induce non-adiabatic changes 
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such as radiation emissions. A mathematically complete perturbation theory 
can obtain not only evolution equations for parameters of solitons analytically 
but also a formula for calculating the perturbation-induced radiation emission. 



With vanishing boundary conditions, — > as |x| — > oo, the derivative non- 
linear Schrodinger(DNLS) equation, 



is an integrable model describing small amplitude nonlinear Alfven waves in a 
low- (3 (the ratio of kinetic to magnetic pressure) plasma, propagating strictly 
parallel to the ambient magnetic field. Here u, x, and t represent the transverse 
complex magnetic field, the normalized longitudinal coordinate, and the nor- 
malized time, respectively and the subscripts denote partial derivatives [4, 5, 6]. 

The modified nonhnear Schrodinger (MNLS) equation. 



is an integrable model in describing propagation of femtosecond pulses in sin- 
gle mode fibers [7,8,9,10]. For picosecond pulses, it is well known that the 
nonlinear Schrodinger (NLS, the case when s = 0) equation is a good model. 
The MNLS model includes the nonlinear dispersion term(the third term on the 
left) which is one of the several higher order effects becoming more significant 
in the femtosecond region[ll,12]. Here u, z, t, and s denote the normalized 
electric field, the normalized distance, the normalized time measured in the 
frame of reference moving with the pulse, and the relative amplitude of the 
nonlinear dispersion, respectively. We assume s > in this letter because the 
case for s < can be obtained by a simple transformation t — > —t. The DNLS 
equation and the MNLS equation are connected by a gauge-like transforma- 
tion(sce, e.g. [13]). The soliton solution for the DNLS with vanishing boundary 
conditions was found by the inverse scattering transform technique [14] and the 
soliton solution for the MNLS can be found by the gauge-like transformation. 

Perturbation theory for the DNLS/MNLS solitons was first developed by a 
method based on the Riemann-Hilbert problem[15], which had minor errors 
in evolution equations for the initial center and phase and gave no correction 
beyond the adiabatic approximation. Later a direct perturbation method[16] 
corrected these errors and obtained a formula for calculating the perturbation- 
induced radiation emission. With well established perturbation theories, in- 
fluence of higher order effects and methods of controlling soliton shapes and 
frequency shift against perturbations can be studied. However, evolution equa- 
tions for the perturbed MNLS solitons are much more complicated than those 
for the perturbed NLS sohtons[l,2]. Problems solved by perturbation theories 



iut + Uxx + i{\u\^u)x = 0, 



(1) 



iQz + ^Qtt + is{\q\'^q)t + \q\'^q = 0, 



(2) 
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for DNLS/MNLS solitons[15,16] were limited to some relatively simple cases 
so far[15.16.17]. Symbolic technique seems helpful for further studies. In this 
letter, we develop a symbolic technique, which is effective for usual perturba- 
tions, to calculate adiabatic evolution equations for parameters of perturbed 
DNLS/MNLS solitons automatically. In the MNLS model, only one higher 
order effect beyond the NLS model, the nonlinear dispersion, is considered. 
Other higher order effects, such as the intrapulse Raman scattering and the 
third-order dispersion, should be treated as perturbations. Evolution equa- 
tions for parameters of the MNLS soliton under some control method, such 
as narrow-banded filter and nonlinear gain, will be valuable in estimating 
these methods. In section 2, we rewrite those evolution equations obtained 
in Ref.[16] in slightly simpler forms. In section 3, we describe the symbolic 
technique. In section 4, we calculate evolution equations of the MNLS solitons 
under intrapulse Raman scattering with local approximation, the third-order 
group velocity(GVD) dispersion, and narrow band filters with nonlinear gain. 
In small s limit, all these results approach those for the NLS solitons in the 
literature, as generally shown in Ref.[16]. 



2 Adiabatic evolution of pEirameters of the perturbed DNLS /MNLS 
solitons 

2.1 Perturbed DNLS solitons 

The DNLS one-sohton[14] can be rewritten as 



Us{x,t) 



2Asin(2/5) 



cosh(^ - ip) 



(3) 



cosh^ie + ip) 



where 



e{x,t) 



2r){x + 4:^t-Xo), 



(4) 



^{x,t) 



(5) 



^ = A^cos(2/5), = A2sin(2/3), < /? < 7r/2. 



(6) 



Its amplitude is 



A = 4Asin/3. 



(7) 
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Simple calculations yield its energy, 



+00 

J dxlusl"^ = 8/3, (8) 



-00 



and the FWHM width, 
1 

A2sin(2/3) 



w = . . ^.^g. In (cos/? + ^l + cos2/?) . (9) 



In presence of perturbations, the zero on the left hand side of Eq. (1) should 
be replaced with the perturbation function ir{u). As /3 is proportional to the 
energy, considering evolution of /3 and A, instead of ^ and 77 in Ref. [16], is 
more convenient. Evolution equations for parameters of the DNLS solitons in 
Ref. ([16]) are reformulated as. 



df3 1 



00 



dt 16r] 



deR+{e,t), (10) 



dA A 



dt I677 



J de tanh{e + ip)R+{e,t), (11) 



^ = #3 / d(^o'-^^-^R.ie,t) - / deR.ie,t), (12) 

dt 16r?3 J cosh(9 + i0) ^ ' ^ 32r?2 J ^ ' ^' ^ ^ 



16773 J cosh{9 + ip) ' ' ' 32?7- 



= 4 r cosh(^-.^) 



8773 i ' ' ' I6772 J cosh{9 + ii3) 



where 



R±{e,t) ^ Uo{e)[ro{e,t) ±ro{-e,t)], (14) 
Uo{9) = Usexp{iip), ro{9,t) = rexp{iip), (15) 

and the bar stands for complex conjugate in this letter. 
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2.2 Perturbed MNLS solitons 



The MNLS one-soliton is 

q,{t,z) = -2Asin(2/3)^^^exp(-i<^), (16) 



where 

9{t,z) ^4:sr]{t-to- z/v), (17) 

^^t, z) = - ^s\e + rf)z - i + ^ + (^0, (18) 

^ = A2cos(2/3), = A2sin(2/3), < /3 < 7r/2, (19) 

and V = s/(l — 4s^^) is its velocity. Similar to the DNLS soliton, we have 
amplitude, 

A = 4Asin/3, (20) 



energy, 

+00 



J dt\q,\'^^, (21) 



and the FWHM width, 
1 

2sA2 sin(2/3) 



r = - •"^(2/3) ^'' (cos/?+^l + cos2/?) , (22) 



of the MNLS soliton. Also, in presence of perturbations, the zero on the right 
hand side of Eq. (2) should be replaced with the perturbation function ir{q). 
Evolution equations for MNLS soliton parameters in Ref. [16] are reformulated 
as. 



dA . A 



hoo 



dz I677 



i— J detanh{e + i(3)R+{e,z), (24) 



5 



^ = ^ T^gg '"t^f/^T? i^-(^, z) 7d9R.{e, z), (25) 

dz 325773 J cosh(e + ip) ^ ' ^ 64^772 J V > ^ V y 

—00 — c» 

— cx) —00 

where 

R±{0, z) = ^[ro(^, ± ro{-e,z)l (27) 
?o(^) = ?sexp(i(^), ro(6', 2;) = rexp(i(^). (28) 



3 A symbolic technique to calculate integrals in evolution equa- 
tions for soliton pctrameters 



In general, integrals in evolution equations for perturbed DNLS/MNLS soli- 
tons are rather complicated. Even modern symbolic softwares can not always 
tackle them directly. However, taking the perturbed MNLS solitons as an ex- 
ample, most perturbation functions can be expressed as 

n 

r{q)^Y.^u{\qmq, (29) 

k=0 



where r^dgp) are complex functions[ll,12]. In what follows we will show that 
for this category of perturbation functions integrals in the evolution equations 
can be systematically solved by the technique of residue theorem. 

For perturbation functions in category of Eq. (29), there are only four types 
of integrands in the evolution equations needed to be tackled. Having been 
continuated to the whole complex plane of 9, they are 

(1) f{e) withf{e + t7r) = -f{e), 

(2) g{e) with g{e + ITT) =g{e), 

(3) F{9)=9f{e), 

(4) Gie)^egi0). 

Within the closed path shown in Fig.l, all of them have and only have a pair 
of singularities 9± = i{7i/2 ± (3). Integrations of them on the two vertical line 
segments are obviously zero. Using the residue theorem on the closed path 
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shown in Fig. 1, we get 

oo 

J f{e)de = m{Res[/(^+)] + Res[/(^_)]}, (30) 

— oo 

and 

oo oo 

j F{e)de = j f{e)de + m{Res[F{e+)] + Res[F(^_)]}. (31) 

— oo — oo 

For g{6) and G{6), by introducing two auxiliary functions, 

M^,p)=exp(^p%(^), (32) 

and 

H(e,p) ^exp{ipe)G(e), (33) 

in which p > 0, we also have 

J ^ ' ^0+ l-exp(-p5r) ^ ' 



OO ^ 

J G{e)de= lim j-TT^ 



Res[/i(^+, p)] + Res[/i(^-, p)] 
2sinh^(p7r/2) 



+ ^.,^ Res[//(gH-,p)] + Res[//(g-,p)] 1 _ 
1 — exp(— p7r) J 

Therefore, despite the fact that these integrals may not be tackled directly, 
calculation of them comes down to calculation of residues and limitations 
which can always be done symbolically by modern commercial mathematical 
softwares. 
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4 Examples of perturbed MNLS solitons 



4-1 Intrapulse Raman scattering 



In local approximation, the intrapulse Raman scattering was described by a 
perturbation function [11, 12], 

r{q) = -(TRq{\qW (36) 



where is a constant. We get -R+ = 2qQro, R- — 0. This problem is so 
simple that we don't need the symbolic technique in the preceding section. 
Direct integrations yields 

^ = -^saRA^(2 + cos^ 2(3 - Q(3 cot 2(5) (39) 
dz 3 



It is obvious that d^/dz monotonically decreases for all /? e (0, 7r/2). The 
MNLS soliton perturbed by the intrapulse Raman scattering retain some sim- 
ilarities with the NLS soliton. They both have no shift in initial position and 
phase. Their energies are not perturbed. They all have frequency redshifts(but 
in different rates). However, the MNLS sohton has a decrease in amplitude A 
and an increase in width r while keeping its energy unperturbed. 

4-2 Third- order GVD dispersion 

The third-order GVD dispersion[ll,12] is described by 

r{q) = adlq, (40) 



where a represents its strength. We have i?+ = 0, i?_ = 2goro, 



and, using the technique developed in the preceding section. 
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dtp 
dz 



a[3s-2-12(/37/ + 4C-/3eV') 



(42) 



d(po 



(43) 



Within adiabatic approximation, influences of the third-order GVD dispersion 
on the MNLS sohton is similar to those on the NLS soliton: the amplitude, 
width, and the main velocity v are unchanged while inducing shifts on 
and (fQ. Numerical simulation showed that the third-order GVD dispersion 
induces radiation emission from the NLS soliton [18] and the problem needs 
to be solved by a perturbation theory beyond all order [19]. With a similar 
simulation, one can find similar radiation emission from the MNLS soliton. 
This means that the problem may also need a perturbation theory beyond all 
order. 

4 .3 Narrow hand filters with nonlinear gain 

For NLS solitons, periodic insertion of narrow band filters was shown to be 
effective in reducing the frequency shift of solitons. But it was found that 
filters may induce background instabilities. M. Matsumoto et. al. suggested to 
utilize nonlinear gains to suppress such instabilities[20]. For further studies, 
we calculate the adiabatic evolution equations of a MNLS sohton controlled 
by both of filters and nonlinear gains here. The corresponding perturbation 
function is 



where F is an excess gain, /t > represents the strength of the filter, pi and 
P2 are nonlinear gain coefficients. We get = 2goro, R- = 0, and, with the 
symbolic technique in the preceding section. 



r{q) = Vq + k— + pi\qfq + P2\qfq, 



(44) 




dipo 



dz 



(45) 



dp 
dz 



2s-'^k(3 - ?>2s'^K{^(3rf - ArjC + 9/3^^) 
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+16p2(2/3772-3?7e + 6/3f)l, 



(46) 



dA 
dz 



+ 



77 



+12pi(-377e + 6/3e' + 2/3?7') 
+ 32p2(477^ - 18/377^^ + 1577^' - ^m^) 



(47) 



5 Summary and discussion 



In this letter, we develop a symbolic technique to calculate the adiabatic evolu- 
tion of perturbed MNLS solitons. Evolution equations under intrapulse Raman 
scattering, third-order GVD dispersion and narrow band filters with nonlinear 
gains are calculated with the technique. As s ^ 0, to keep parameters in the 
MNLS soliton physically meaningful, ^ and 77 must be 

^ A* 1 u 



The MNLS soliton approach the NLS soliton, 

gs(i, 2;) — >• — 2i/sech6'exp(— i<^), (49) 

where 

e = 2u{t-i), i=-2iiz + to, (50) 
if = 4:ii{t -i) + (p, (p = -4(/i^ + iy'^)z + v3o- (51) 



In Ref.[16], it was generally shown that evolution equations for perturbed 
MNLS solitons approach those for perturbed NLS solitons in small s limits. 
Expanding results in the preceding section near s = 0, one can find they do 
approach their corresponding results for NLS soliton in the literature[ll,12,20]. 
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Fig. 1. The path for integrals in section 3. The crosses are the singularities. 
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